Enhanced subthreshold e + e production in short laser pulses 
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The emission of e + e~ pairs off a probe photon propagating through a polarized short-pulsed 
electromagnetic (e.g. laser) wave field is analyzed. A significant increase of the total cross section of 
pair production in the subthreshold region is found for decreasing laser pulse duration even in case 
of moderate laser pulse intensities. 
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The history of the study of e + e~ production in Y7 
interaction starts with the pioneering work by Breit and 
Wheeler [l[ published in 1934. About thirty years later, 
Reiss and Narozhnyi, Nikishov and Ritus [H, HJ have 
analyzed the e + e~ emission off a photon 7' propagating 
in the field of an intensive polarized monochromatic elec- 
tromagnetic (cm.) plane. The e + e~ production prob- 
abilities were found using the non-perturbative Volkov 
solutions for the electron and positron wave functions Q . 

If one identifies the external e.m. field with a laser 
pulse then most of the early work considers long last- 
ing pulses where the temporal shape can be neglected. 
We denote this approach as the infinite pulse approxi- 
mation (IPA). In IPA, electrons and positrons e + be- 
come quasi-particles with effective quasi-momenta and 
effective (dressed) masses. Differential and total proba- 
bilities of the e + e~ pair emission depend on the reduced 

strength of the e.m. field A», f 2 = -^f 1 = 
where M e is the electron mass (we use c = h = 1, 
e 2 /47r = a = 1/137). Furthermore, the dimensionless 
variable £ = ^f 1 - is introduced, where s is the square of 
the total energy in the center of mass system (c.m.s.) 
of the Breit- Wheeler process 7' + 7 —> e + + e~ and 
Sthr = 4M 2 is its threshold value. The Ritus variable 
is then defined by K = 2£/f Q. The case of £ > 1 cor- 
responds essentially to multi-photon processes. Within 
IPA, the minimum number of photons 7 in the reaction 
j'+nj —> e + +e~ is defined as n min = 7(C) + 1, where 
is the integer part of £. First evidence of the multi-photon 
Breit- Wheeler process with £ = 3.83 and 0.1 < £ < 0.35 
was detected at SLAC in the E-144 experiment [y], where 
the application of IPA is justified since the used laser 
pulses contain around 10 3 cycles in a shot. 

The rapidly evolving laser technology Q can provide 
the laser power up to 10 24 - 10 25 W/cm 2 in near future 
which is sufficient for the formation of positrons from cas- 
cade processes in the photon-electron-positron plasma 13- 
10 1 generated by photon- laser lldHl3| . electron- laser [1J, 
15 1 or laser-laser interactions [l6l \l7\ (see [lH for sur- 
veys). The next generation of optical laser beams are 
expected to be essentially short (femtosecond duration) 
with only a few oscillation of the e.m. field in the pulse 



to be expected at ELI and CLF [Io[ facilities. This 
requires the generalization of the IPA multi-photon pro- 
cess 7' + 717 — > e + + e~ to a finite pulse duration. For- 
mally, this generalization may be done in a straightfor- 
ward manner by substituting the expansion in Fourier 
series into Fourier integrals with taking into account the 
Volkov solution for the finite wave field. In practice, an 
evaluation of the total cross section requires the calcula- 
tion of five-dimensional integrals with rapidly oscillating 
integrands which is rather demanding. Therefore, pre- 
vious considerations are often restricted to the analysis 
of the three-dimensional differential cross sections, see 
for example (l2j for finite beam size effects in e + e~ pair 



production (cf. also 2l[ and references therein). 

The aim of the present Letter is to elaborate a method 
for the calculation of the total cross section in the sub- 
threshold (multi-photon) region accounting for the effect 
of finite laser pulse duration in e + e~ pair production off 
a probe photon. We denote such a process with a finite 
pulse and plane wave fronts as finite pulse approximation 
(FPA). In this case, the in/out fermion states refer to the 
vacuum. Moreover, due to the modulation of the pulse 
envelope function, the power spectrum contains frequen- 
cies > uj (see below) which enhance the pair production 
in the subthreshold region even for moderately strong 
laser intensities. 

We consider the e.m. four-potential A ~ (0, A) in FPA, 
depending solely on the invariant phase (f> = k ■ x, 



A(0) = /(</>) (ai cos</> + a 2 sine 



(1) 



where |ai| = | a.2 1 = a, aia 2 = for circular polariza- 
tion. We employ here the envelope function /(</>) = 
1/ cosh(</>/A), where A = tt^ = irN, and N charac- 
terizes the number of cycles in a pulse; tq = 2it/uj is the 
time of one cycle for the laser frequency u>. Thus, r is 
the time scale of the pulse duration. The case of pulses 
obeying lot 3> 1 has been analyzed in p2j . 

Utilizing the e.m. potential (JXJ) in the Volkov solutions 
leads to two significant modifications of the transition 
amplitude. Besides physical asymptotic momenta and 
masses, the finite time r requires Fourier integrals in the 
integrand of invariant amplitudes, and the discrete har- 
monics become continuous. Thus, the S matrix element 



2 



of the process 7' — s- e + (j) + e (7), where e ± (7) refers to 
Volkov states in the field (U}, is expressed as 



where the transition matrix M(l), similarly to the case of 
the non-linear Compton effect |23l - [2rj as a crossed chan- 
nel of the pair production, consists of four terms 



M(l) = M(m) C {m) (l) 



(3) 



m=0 



where 



C( m )(0 = — / d0 x (m) (</>) e a *- iVW . (4) 

Here, X (m) = (1, / 2 (0), /(0) cos 0, /(</>) sin 0) with m = 
0,1,2,3 and 



V{4>) = zV Q {4>,<Po)-i%u I ^'/ 2 (0'), (5) 

r</> 



M4>, M 



#' cos(0'-0 o )/(<//), 



(6) 



where w = (fc-fc') 2 /(4(fc-p)(/cy )), z = 2^ A /u(u - u)/"o, 
Uo = l/C- The angle O is related to the azimuthal 
angle of the positron in the e + e~ rest frame by 0o = 
4> p + 7r and can be determined through invariants <X\^ = 
e((cti,2 -p) ~ (ai,2 -p')/(k -p')) as cos0 o = a-i/z, 

sin</>o = a.%jz. Here, p e - = p' ~ (pq, p') and p e + = p ~ 
(poiP)- The transition operators M^ 2 ' 3 ' are the same as 
in IPA [H, while the operators M (04) = u p > Ai^ ' 1 ' v p , 
read now 



4(fc-p)(fc-p') : 



(7) 



where u„i and u p are the free-field Dirac spinors of the 
outgoing electron and positron, respectively; e' is the po- 
larization four- vector of the probe photon 7' with four- 
momentum k' ~ (o/,k'), and fc ~ (^,k) is the four- 
momentum of the e.m. (laser) field <[TJ) . Fcynman's slash 
notation is employed, e.g. J/l = A • 7, as four-product 
with the Dirac 7 matrices. The integrand of the function 
C(°> does not contain the envelope function and needs a 
rcgularization, e.g. using a prescription given in Ref. [23j 



Jlc/>-iV(4>) 



X 2: cos 



0o)/W>)-eW 2 (0)) . (8) 



The probability is normalized to some time unit. In 
IPA, one can use the time of one cycle, tq. In FPA, a 
proper time unit is provided by the pulse width, which is 
N times greater, r = Ntq, where ./V is the number of the 
cycles in a pulse. Therefore, for a convenient comparison 



of IFA and FPA results, the latter one is scaled by 1/N. 
Thus, the probability of the e + e _ pair emission reads 



W 



FPA 



2tt 



dl 



iio du ^mdd, ( 9) 



w(l^,u,cj )p ) = (2u a + l)\C^(l)\ 2 

+e 2 (2 U -i)f|c (2) (0l 2 + |c (3) (0l ; 



(10) 



+ReC<°>(J) U 2 C«(0- ^[a 1 cW(l) + a 2 C&(l)} 

with uq = l/C,. This expression will be used below for 
direct numerical evaluations of the probability. 

Inspection of the functions V{4>) and C^ m \l) shows 
however that Eq. (fTO]) may be simplified to get, in some 
cases, a more suitable analytical expression for w(l). In- 
tegrating by parts, the function Vq(4>,4>o) might be ex- 
pressed in the following form 



W, 0o) = sin(0 - O ) + 0(A), 



(11) 



where 0(A) = sin(<// - ^ )/ / (^ / )#' is a rather 

small contribution for a finite pulse duration A = ttN 
with N > 2 because of (i) the factor 1/A and (ii) the 
derivative f'(<j>) in the integrand has a maximum value 
at the boundaries of the pulse with 4> ~ 0.9A, where this 
function is suppressed. In fact, the numerical evaluation 
shows that the contribution of 0(A) can be omitted (we 
find 1 0(A) I < 0.1 (0.05) for A = 2tt (5tt)). This approx- 
imation allows to express the basic functions C^ m '(l) via 
new functions Yi and Xi 

Crt°>(l) = YK^e'^ , C {1) {1) =X l (z)e a *° , (12) 
C [2 \l) = ^(ywWe^ + li-iWe" 1 -^) , 
with 



TO 
Yi(z) 



1 

2^ 
1 

2^ 



/(2) (V> + <fo) e a ^- lz f^+<t> a ) sin * , 



/<"%) = /™(0)expK 2 C«tanh^] 



(13) 



The function C( 3 ) emerges from C^ 2 ' by the substitutions 
1/2 — > l/2i and sign " + " between two terms in the 
bracket to "— ". In the last line, "to" (=1, 2) is a label on 
the l.h.s., while on the r.h.s. it is the power of the envelope 
function, as follows from Eqs. (j4j, ([SJ; the exponential 
term results from an analytic evaluation of the last term 
in Eq. (5) for the chosen envelope function. 
The partial probability w(l) in Eq. ((TU)) reads 

w(l,^u, ( p p ) = 2Y l 2 (z)+e(2u-l) 

x (y^z) + Yu-i(z) - 2ReY l (z)Xt(z)) (14) 



3 



which resembles the expression for the probabilities w n 
in case of IPA (cf. Ref. Q) arising upon the substitu- 
tions Jdlw(l) -> ^ n w n , Y? -> J 2 , y ; 2 ±1 -> J 2 ±1 , 
ReYi(z)Xi(z) J 2 with Bessel functions J„. 

In the case of small field intensity, £ <C 1, implying 
z < 1, and denoting / = n + e, where n is the integer 
part of I, one can use the following decomposition 



Y, 



1 

2^ 



etye , ' , *- fa/w+ * o) ™^/(V> + 0o) (15) 



1 



and analog for the function Xi(z) with the substitution 
jk+i _^ / fc + 2 . The dominant contribution to the integral 
with a rapidly oscillating integrand stems from the term 
with k = n. which results in 



Yi 



k+<L 



X, 



k+e 



2 k k\ 
z k 
2 k kl 



(16) 



where the function fp .(e) is the Fourier transform of the 
function f k {ip). For the above envelope function it can 
be calculated analytically using the theory of residues. 
Results of the leading orders n = 0, 1 are 



Y 0+e (z) 

Yl+e(z) 



A -vrlelA/2 

e - ie </> 

1 + c - 7r l £ l A 

z A 2| £ | e -^| e |A/2 

2 i^T^RA 
zA(A¥ + l)c- 
4~ 



r|e|A/2 



1 



7r|e|A 



(17) 



The representation of Eq. (|16|) evidences (i) a fast de- 
crease of Y 2 +e with increasing |e| and (ii) the <j> p depen- 
dence disappears in Y 2 +e and X% +e . This allows to ex- 
press the integral over dl in ([9]) in a form useful for a 
qualitative analysis: 



( Wn\ 



w 



FPA 



oo 



E 

n— TT-o + 1 



ie r = no,e,Z,u) 

C-n h U 3 ' 2 Vu - 1 

\e Hdu^M (18) 



U 3 ' 2 y/u — 1 



with uo = (n + e)/Cl n o = 1 f° r C ^ 1; an d n = /(£) for 
C > 1; The lower limit in integral over de in the second 
term reads v = Q—n for £ > 1 and n = no + 1, and f = — 1 
in other cases. This equation shows that, contrary to IPA 
where at given £ > 1 (i.e. below threshold, s < sthr = 
4M 2 ) only harmonics with n > /(£ + 1) contribute, in 
FPA the harmonic with n = /(£) also contributes. 

Consider, as a check of the normalization, the pair pro- 
duction above threshold with £ = 1 — Ss e /s < 1, where 
5s e is the energy excess 5s e = s — Sthr- Utilizing the ex- 
plicit expressions (|17[) for the leading contribution Yi+ e 



in (jl4[) one can get a relation between emission probabil- 
ities in IPA (cf. H) and FPA: 



W 



FPA 



IF IPA (?i = 1,£, ui)Z(A,£), (19) 

e -7r|e|A 

de— 

><- 



A 2 f 1 
I(A,0 = — / de— 

{ 7 n y c _i 1 



a -7r|e|A^2 



(20) 



where u n is an effective value of u in nth term of Eq. (|18[) . 
The dependence of VI 7 on ti„ is rather weak compared to 
the dependence on £ and can be disregarded. Thus, in 
the limit irA6s e /s S> 1, IPA and FPA practically coincide 
since Z(A, £) ~ A[nN(l + e-^-W^j-i ~ ^ = 1. 

Consider now the case of subthreshold pair production 
with £ = 1 + <5s;/s > 1, where 5s; = s t i ir — s is the "lack 
of energy" . The probability has the following form 



W 



FPA 



W IPA (n = 1) 



CW lpA (n 



2) 



(21) 



-nA5 Sl /s^ 1 + e -nA6 Sl / S ^ and Q 



with Xi(A,5si/s) ~ e 

(1/-7T 2 ) j*^Q_ 2 ) x 2 exp(— x)(l — exp(— a;))~ 2 d;r ~ 2/3 for 
<5s//s < 1 — 0.65/iV. The terms in r.h.s. of (f2~Tj) are meant 
to have the same functional dependence on £ and u 12 
as in IPA. One can expect a significant enhancement of 
pair production for the short pulse because the proba- 
bility of single photon events (n = 1) is much greater 
than the probability of the two-photon events (n = 2): 
W lFA {n = l)/W lpA (n = 2) ~ £~ 2 > 1. When the 
length of the pulse increases the contribution of the first 
term in Eq. ([21]) decreases exponentially due to X\, and 
the prediction of FPA approaches to the IPA one. 

The probability and the cross section are related to 
each other ^ as dW = 2[loM It? j {^u)\do . The to- 
tal cross section of e + e _ production is calculated using 
Eqs. ([9]), (fT4]) and (fl6|) . The cross sections are exhib- 
ited in Fig. [1] as a function of yfs in the threshold region 
for finite pulses with A = irN. The left and right pan- 
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FIG. 1: The total cross section of e + e~ pair production as a 
function of the total energy in c.m.s., y/s, for the finite pulse. 
Notation is given in the text. 

els correspond to £ = 0.01 and 0.1, respectively The 
dashed and thick solid curves are for N = 2 and 5, re- 
spectively. The thin solid curve is the IPA result. The 
thin dashed curve, labelled by "B-W" , corresponds to the 
Breit- Wheeler process [l| practically coinciding with the 
lowest harmonic (n = 1). One can see that in the sub- 
threshold region, s/s — 0.85 - 1.02 MeV, the cross section 
for short pulses is significantly greater than in IPA and 
the difference may reach one or two orders of magnitude 



4 



for £ = 0.1 and £ = 0.01, respectively. When £ and/or 
( increase, the contribution of higher terms with n > 1 
becomes finite that brings an additional (increasing) de- 
pendence on A (cf. Eq. (jTGj) . 




FIG. 2: Left panel: As Fig. [T] (right panel) but for a wider 
region of y/s. Right panel: The total cross section of the 
e + e~-pair production as a function of £ 2 at three values of 
Imm = C — 0-5, 1-1 and 3.8. Notations as in Fig. [T] 

The total cross section in a wider region of y/s is exhib- 
ited in Fig. [2J left panel. At x/s ~ 0.55 MeV the multi- 
photon events with I > 4 become important. In general, 
the total cross section in FPA has also the step- like struc- 
tures similar to IPA. However, a decrease of the pulse 
duration leads to a smoothing. One can also see some 
enhancement of the cross section for a short pulse with 
N = 2 compared to the case of a longer pulse with N = 5. 
The total cross sections of the e + e~ pair production as 
a function of £ 2 at three values of l mm = C = Sthr/s are 
presented in Fig. [2J right panel. The case of £ = 0.5 cor- 
responds to the production above the threshold. Here, 



the predictions for IPA and FPA coincide. Examples of 
C = 1-1 and 3.8 correspond to the subthreshold produc- 
tion. In the first case, we are slightly below the threshold 
and one can see a large difference between predictions 
for pulses with N = 2 and 5, which has been explained 
above. The last example (£ = 3.83) corresponds to the 
kinematics of the SLAC E-144 experiment. In this case 
the predictions of IPA and FPA arc qualitatively simi- 
lar with some enhancement for a shorter pulse duration. 
Finally note that we do not take into account radiation 
reaction effect discussed in [|| and recently in Ref. [28| 
because it influences the fermions in the final state and is 
not expected to change significantly the total e + e~ yield. 

In summary, we have considered the total cross section 
of e + e~ production off a probe photon interacting with 
a semi-intensive short laser pulse in the subthreshold re- 
gion defined by multi-photon interactions. We find a non- 
trivial dependence of the cross section (production prob- 
ability) on the pulse duration. Just below the threshold 
of the weak-field Breit- Wheeler process, the short laser 
pulses increase the cross section up to two orders of mag- 
nitude relative to a monochromatic plane wave. This 
effect must be taken into account in the evaluation of 
e + e~ pair production in cascade processes produced by 
high-power laser fields. 
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for the kind hospitality. 
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